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SYSTEM AND METHOD FOR 
ANALYZING A PATTERN IN A TIME-STAMPED EVENT SEQUENCE 



FIELD 



[0001] The present disclosure relates to data processing, and more 
particularly, to the field of data mining of time-stamped events in a temporal 
sequence and identification of a surprising pattern in the sequence. 



10 



BACKGROUND 



[0002] With prevalent application of computer technology to 
transactions and business operations, a very large amount of operational event 
data is being generated and stored in databases. In many applications, the 
stored event data includes a time-stamp that provides for the identification of the 
1 5 time of occurrence of the event. While a very large amount of sequential data is 
stored in databases, generally only a limited amount is analyzed due to the high 
cost in mining the data. As such, many patterns of interesting events remain 
hidden. 



20 define relationships between two or more data events. One such method for 

identifying patterns of data events within sequences is data mining. Data mining 
generally is the extraction of knowledge or patterns from data in databases or 
other information repositories. In contrast to simple database searches, data 
mining finds each sequence with at least one pattern satisfying the constraint. 

25 [0004] A large volume of event data has been collected in applications 

such as maintenance records and web click applications. For example, a 
sequence of maintenance events may include a list of operational or failure 
events ordered by time of occurrence. Each event data item may include an 
identification of the particular event, the type or categorization of the event, and 

30 the time of occurrence. While events may be ordered in time, their contents have 
no ordering and are not easily compared to identify a trend or pattern. 

[0005] One example of maintenance events stored in a temporal 
sequence is operational events associated with an aircraft. Each sequence may 
be associated with a particular aircraft and two or more sequences may be 



[0003] It is generally desirable to identify patterns of data events that 
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associated with a fleet of aircraft. Such aircraft maintenance eventsequences 
are different than a simple time series. Maintenance events occur irregularly. 
Often, some time periods do not contain an event and other time periods 
containing two or more "overlapping" events. 
5 [0006] For sequences containing temporal maintenance events, data 

mining methods and systems are generally designed to identify events that 
precede a hardware failure or maintenance event. The identified events usually 
are within a monitoring time range at which intervention may be possible to 
prevent a failure or to reduce a cost associated with the next event. Such 

1 0 patterns of events are generally subject to ordering and temporal constraints. 
Pattern occurrence is a fundamental concept in the sequential pattern data 
mining problems. As such, data mining methods may identify a pattern that 
forecasts a target event such as a failure or operational event. 

[0007] A data mining task may include discovering sequential patterns 

15 among events, i.e., co-occurrences of multiple events and some ordinal or 
temporal relationship among them. The discovered patterns may be then 
interpreted as rules. An example pattern is an Engine Oil event followed by an 
Automatic Flight event in one to three days. This pattern can be interpreted as a 
sequential rule such as: If an Engine Oil event occurs within one to three days, 

20 an Automatic Flight event will occur. 

[0008] As such, sequential pattern mining methods generally utilize 
pattern occurrence identification techniques. Pattern occurrence identification 
methods have been developed for sequential pattern discovery, filtering and 
ranking. In these methods, generally the recurrence of a pattern within the same 

25 sequence is ignored. Additionally, frequent pattern mining generally considers 
multiple sequences and often ignores pattern recurrences within a single 
sequence. However, the number of pattern occurrences in a single sequence 
can provide valuable insight especially for applications having long sequences 
each containing many events such as maintenance record data. For example, 

30 airplane maintenance records are usually kept for the life time of each airplane, 
and many patterns naturally repeat in the maintenance history of the same 
airplane. The number of occurrences within sequences might indicate problems 



2 



Substitute Specification (Clean Version) 

Boeing Ref. 02-1381 
HDP Ref. 7784-000707 

of a particular airplane (while the number of occurrences across sequences might 
indicate problems of a group of airplanes). 

[0009] One data mining method is a constraint-based mining method 
that utilizes user defined constraints defining the pattern to be mined and 
5 includes classification and association constraints. Another method includes 
distance-based association rules such as the density or number of events in an 
interval and/or the closeness of events in the interval. 

[0010] Another method provides for the discovery of sequences of 
maximum length with support above a given threshold. A sequence is defined as 

10 an ordered list of elements where an element is defined as a set of items 
appearing together in a transaction. This method identified two data mining 
metrics, support and confidence. Support is defined as the extent to which the 
data is either positively or negatively relevant to the rule. Confidence is defined 
as the extent to which, within those that are relevant, the proposal is upheld. 

1 5 [001 1] Another method uses a sliding window on the input sequence to 

obtain a set of overlapping subsequences, and reports the number of 
subsequences in which the pattern occurs. Recurrences within a subsequence 
are ignored. Different numbers of occurrences for a pattern are a function of the 
selected window size. When the window size is large enough, all legitimate 

20 occurrences are considered. However, the same event instances or event 

pattern occurrences may be counted multiple times in multiple sliding windows 
even though there are only two instances of a particular event. The number of 
pattern occurrences increases as the window size increases. However, this 
method is limited as the choice of window size is critical. In addition, the sliding 

25 window approach is static and not very robust. For example, increasing the 
window size introduces a different number of new occurrences for different 
patterns, and thus changes the order of patterns in terms of the pattern 
occurrence or other derived measures. 

[0012] In another method, only the minimal pattern occurrences are 

30 counted. In such a method, an occurrence is identified as minimum if no other 
occurrence can be found in any proper sub-interval of its time span. Legitimate 
occurrences of the pattern that are not "minimal" are ignored. However, a more 
constrained pattern may have more minimal occurrences. As such, such a 
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method produces an unexpected result due, in part, to the exclusion of some 
legitimate occurrences. 

[0013] Another data mining method includes the identification of an 
interesting pattern where events of an episode occur close in time. An episode is 
5 a conjunction of events bound to a given variable and that satisfies unary and 
binary predicates declared for those variables, e.g., a collection of events 
occurring frequently together or partially ordered collection of events occurring 
together. The method distinguishes between serial and parallel episodes and 
between simple and non-simple episodes, where a simple episode contains only 

10 unary predicates and no binary predicates. In this method, a time window is a 

user defined width of time defining how close the events must occur to each other 
within the episode. A window is a slice of an event sequence. An event 
sequence is a sequence of partially overlapping windows. The user may also 
specify how many windows an episode has to occur to be considered a frequent 

1 5 episode. Episodes that occur frequently within a sequence are determined. 

[0014] In yet another method, a number of disjoint occurrences is 
determined. This method addresses discreet events and their relationship to 
each other, but does not allow for time overlapping events within the sequence. 
As such, this method is not applicable to patterns and sequences containing 

20 maintenance events or web transactions that inherently have time overlapped 
events. 

[0015] Each of these methods is limited in their application and 
effectiveness in determining or identifying a pattern within a sequence of time- 
stamped events or categories. Therefore, the inventor of the present method and 
25 system believes it would be desirable for a method and system to effectively and 
efficiently provide for the identification of a pattern in a sequence of time-stamped 
events. The inventor also believes that it would be desirable for a method to 
provide for the identification of surprising patterns within one or more temporal 
sequences. 



4 



Substitute Specification (Clean Version) 

Boeing Ref. 02-1381 
HDP Ref. 7784-000707 

SUMMARY 

[0016] In one implementation, a method determines distinct 
occurrences of a pattern in one or more sequences of time-stamped event 
instances. The method includes determining a maximum cardinality of disjoint 
5 occurrences of the pattern in the one or more sequences. 

[0017] In another implementation of the method, an expected quantity 
of distinct occurrences of a pattern in a sequence of time-stamped events 
assigned to event categories is determined. The pattern includes a first event 
category and a second event category that is within a time gap of the first event 

10 category. The time gap has a minimum time gap and a maximum time gap, and 
the sequence has a maximum time length. The method includes counting 
instances of the first event in the sequence and counting instances of the second 
event in the sequence. The method also includes determining the expected 
quantity of distinct occurrences of the pattern as a function of the quantity of first 

1 5 event instances, the quantity of second event instances, the maximum time 
length of the sequence, the minimum time gap, and the maximum time gap. 

[0018] In yet another implementation of the method, a surprise pattern 
within a sequence of time-stamped event instances is identified. The method 
includes calculating an expected quantity of distinct occurrences of a pattern in 

20 the sequence. The method also includes determining a maximum cardinality of 
the pattern in the sequence. The method further includes identifying the surprise 
pattern as a function of the estimated quantity of distinct occurrences and the 
maximum cardinality. 

[0019] In still another embodiment, a system determines distinct 

25 occurrences of a pattern in a sequence of time-stamped event instances. The 
system includes means for storing the sequence and means for defining the 
pattern. The system also includes means for determining a maximum cardinality 
of disjoint occurrences of the pattern in the sequence. 

[0020] In another embodiment, a computer readable medium includes 

30 computer executable instructions for determining distinct occurrences of a pattern 
in a sequence of time-stamped event instances. The computer executable 
instructions include means for determining a maximum cardinality of disjoint 
occurrences of the pattern in the sequence. 
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[0021] In still another embodiment, a system estimates an expected 
quantity of distinct occurrences of a pattern in a sequence of time-stamped 
events assigned to event categories. The pattern has a first event category and 
a second event category, with the second event category being within a time gap 
5 of the first event category. The time gap has a minimum time gap and a 

maximum time gap and the sequence has a maximum time length. The system 
includes means for counting instances of the first event in the sequence and 
means for counting instances of the second event in the sequence. The system 
also includes means for determining the expected quantity of distinct occurrences 

1 0 of the pattern as a function of the quantity of first event instances, the quantity of 
second event instances, the maximum time length of the sequence, the minimum 
time gap, and the maximum time gap. 

[0022] In yet another embodiment, a computer readable medium 
includes computer executable instructions for estimating an expected quantity of 

1 5 distinct occurrences of a pattern in a sequence of time-stamped events. The time 
stamped events are assigned to event categories. The pattern has a first event 
category and a second event category, wherein the second event category being 
within a time gap of the first event category. The time gap has a minimum time 
gap and a maximum time gap. The sequence has a maximum time length. The 

20 computer executable instructions include means for counting instances of the first 
event in the sequence and means for counting instances of the second event in 
the sequence. The computer executable instructions also includes means for 
determining the expected quantity of distinct occurrences of the pattern as a 
function of the quantity of first event instances, the quantity of second event 

25 instances, the maximum time length of the sequence, the minimum time gap, and 
the maximum time gap. 

[0023] In another embodiment, a system identifies a surprise pattern 
within a sequence of time-stamped event instances. The system includes means 
for storing the sequence of time-stamped event instances and means for defining 

30 the pattern. The system also includes means for calculating an expected quantity 
of distinct occurrences of a pattern in the sequence. The system further includes 
means for determining a maximum cardinality of the pattern in the sequence. 



6 



Substitute Specification (Clean Version) 

Boeing Ref. 02-1381 
HDP Ref. 7784-000707 

The system also includes means for identifying the surprise pattern as a function 
of the estimated quantity of distinct occurrences and the maximum cardinality. 

[0024] In yet another embodiment, a computer readable medium 
includes computer executable instructions for identifying a surprise pattern within 
5 a sequence of time-stamped event instances. The computer executable 
instructions include means for calculating an expected quantity of distinct 
occurrences of a pattern in the sequence. The computer instructions also include 
means for determining a maximum cardinality of the pattern in the sequence. 
The computer instructions further includes means for identifying the surprise 

1 0 pattern as a function of the estimated quantity of distinct occurrences and the 
maximum cardinality. 

[0025] Some implementations and embodiments of the present 
disclosure provide for improved efficiency and effectiveness for mining patterns in 
sequences of time-stamped events. This provides for lower data mining costs 

15 and increases the opportunity for mining data. Some embodiments also provide 
for improved identification of surprising patterns in one or more sequences. 

[0026] Further aspects and features of the present disclosure will be in 
part apparent and in part pointed out in the detailed description provided 
hereinafter. The features, functions, and advantages can be achieved 

20 independently in various embodiments of the present disclosure or may be 
combined in yet other embodiments. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0027] The present disclosure will become more fully understood from 

the detailed description and the accompanying drawings, wherein: 

[0028] FIG. 1 is a sequence of time-stamped categories of events 

5 according to one embodiment of the present disclosure. 

[0029] FIG. 2 is a flow chart illustrating a method of determining a 

maximum cardinality of a pattern in a sequence according to one implementation 

of the present disclosure. 

[0030] FIG. 3 is a flow chart illustrating a method of determining a 

10 maximum cardinality of a pattern within a sequence according to another 

implementation of the present disclosure. 

[0031] FIG. 4 is a flow chart illustrating a method of identifying event 

instances in occurrences within a sequence according to one implementation of 

the present disclosure. 

1 5 [0032] FIG. 5 is a flow chart illustrating a method of determining a 

maximum cardinality of a pattern within a sequence according to yet another 

implementation of the present disclosure. 

[0033] FIG. 6 is a flow chart illustrating a method estimating the 

expected maximum cardinality of a pattern within a sequence according to one 

20 implementation of the present disclosure. 

[0034] FIG. 7 is a flow chart illustrating a method estimating the 

expected maximum cardinality of a pattern within a sequence according to 

another implementation of the present disclosure. 

[0035] FIG. 8 is a flow chart illustrating a method identifying a 

25 surprising pattern in a sequence according to one implementation of the present 

disclosure. 

[0036] FIG. 9 is a flow chart illustrating a method identifying a 
surprising pattern in a sequence according to another implementation of the 
present disclosure. 

30 [0037] FIG. 10 is a functional block diagram of a system for 

determining maximum cardinality of a pattern in a sequence, expected frequency 
and/or surprise patterns from a time-stamped event sequence according to one 
embodiment of the present disclosure. 
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DETAILED DESCRIPTION 

[0038] The following description of the implementations and 
5 embodiments is merely exemplary in nature and is in no way intended to limit the 
present disclosure, its application, or uses. For purposes of clarity, the same 
reference numbers are used in the drawings to identify similar elements. 

[0039] Before describing a system or method for determining the 
maximum cardinality of a pattern in a sequence, determining the expected 
10 maximum cardinality, and identifying a surprising pattern in a sequence, one or 
more concepts associated with a sequence will be introduced and defined. 
Thereafter, a detailed description of various embodiments will be described. 

[0040] Following is a list of symbols and terms used throughout this 
specification. This listing is intended for illustration purposes and is not intended 
15 to be limiting. 

1 . e - an event instance. 

2. A, B, C, D, and E - Capital Letters are indicative of categories of events. 

3. C(e) - a category of event instance e 

4. b = a factor for determining the mean of the expected maximum 
20 cardinality 

5. c - a maximum cardinality of pattern P in sequence s. 

6. d - a parameter that is a function of alternating adjustment factor over a 
range of indices is used to determine the incremental estimation 
parameter ip. 

25 7. g - a total number of groups in P pattern. 

8. group - a set of events that may contain multiple copies of the same 
event categories with a group window size constraint. 

9. h - a secondary index of event i used to go backward from event i. 

10. i - an index of an event in the whole sequence. 
30 11. j - an index number of the group. 

12. k - a loop count index. 

13. / - a data sequence maximum time range. 

14. m - a total number of events in sequence s. 

15. n - a number of event categories in a pattern, where nj is the number of 
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event categories in the j-th group in a pattern P. 

16. P - a pattern that is a collection of event categories with structural and 
temporal constraints. It is an ordered list of groups with minimum and 
maximum time gap constraints between any two consecutive groups. 
The gaps and window sizes are integers indicating time differences as 
specified in a time unit. Different gaps and window sizes may be present 
in the same pattern. If a group contains a single category, its window 
size is 0 and is omitted with the colon separator. Pattern P may, for 
example, be designated {A}a - (3{B}. 

17. p - a probability of the occurrence of one event or category in relation to 
a second event of category. The probability p is approximately equal to 
the remainder of the maximum time gap minus the minimum time gap 
divided by the length of the sequence, e.g., p=(p - a)/(/J . For any given 
instance of event A, there is a probability p that any event instance of B 
together with the event instance of A are an occurrence of the pattern P. 
For example, for a sequence having a length / of 1 ,000 seconds and a 
pattern P of {A}1-10{B} where B is followed within 1 to 10 seconds of A, 
the probability is equal to 1-q=10/1000=.01=1 percent. As such, for any 
given event or category A instance within the sequence, there is a 1 
percent change for any B instance such that these two instances are an 
occurrence of pattern P. 

18. q - a parameter equal to one minus the probability p. 

19. S - a set of a disjoint pattern. 

20. s - a sequence of temporal events. 

21 . T(e) - a timestamp of the event instance e in sequence s. 

22. x - a number of event instances of category or event A in sequence s [1 , 
/]. 

23. X x , y - a random variable of the number of occurrences of pattern {A}a- 
0{B} having x occurrences of A's and y occurrences of B's and that 
where y > x > 0. 

24. y - a number of event instances of category or event B in sequence s [1 , 
/)• 

25. z - an estimation index. 
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26. aj - a minimum time gap between the j-th group and the (j + 1 )-th group. 

27. ft - a maximum time gap between the j-th group and the (j + 1 )-th group. 

28. y - an estimation adjustment parameter. 

29. 6 - a base estimation of the mean of the expected maximum cardinality. 
5 30. r|x - a number of instances of type C x events. 

31 . A - an alternating adjustment factor. 

32. p - a bound to the mean of the expected maximum cardinality where p " 
is the minimum bound and p + is the maximum bound. The real mean is 
denoted a p r 

10 33. Sj )j>b - a probability that Xjj equals b, i.e., ^ Jib := Pr(Xjj = b). 

34. p - an estimation precision objective that is a fraction of real mean u r 
such that the estimated bounds are tighter than this fraction, i.e., the 
upper bound minus the lower bound must be less than p. 

35. a - a loop counter where a is one plus the number of loops or 
1 5 recurrences of the method. 

36. q> - an estimation coefficient. 

37. - an incremental estimation parameter. 

38. coj - a window size of the j-th group. 

20 [0041] An event instance e has a unique ID, a category name, and a 

timestamp. As indicated by sequence 100 in Fig. 1, timeline 102 indicates time 
periods from 0 to 26. Sequence 100 is a set of event instances of categories are 
indicated as A, B, and C. For example, event instance category C 104 occurred 
at time 0. The next event category A instance 106 is time stamped at time 1 1 . A 

25 second event category A instance 108 is time stamped at time 13. Event 

category B instances occur at time periods 14, 15, and 16 as indicated by 110, 
112, and 1 14. A third event category A instance 116 also occurs at time period 
16 along with category B instance 1 14. A second category C instance 118 
occurs at time period 25. Each of event category A, B, and C as indicated by 104 

30 to 1 18 comprise sequence 100. As noted, in a sequence multiple event 

instances with the same category for A, B, and C may occur. Additionally, the 
sequence may include multiple instances of the same or of different categories at 
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the same timestamp such as instance B 1 14 and instance A 1 16 at time period 
16. 

[0042] The method of the present disclosure may be used to analyze 
sequence 100 or other time stamped events in any sequence of time stamped 
5 events or categories. Sequences with different time units and difference time 
origins may also be analyzed with this method and system. Integer numbers are 
used in Fig. 1 to denote the event timestamp T, however, any numerical counting 
of a timestamp may apply. 

[0043] A sequential pattern P may be described by one or more 
10 characteristics that may include those identified in Table 1 . 



Term 


Definition 


(pattern) 


(group)((min gap) "-" (max gap) 
(group) )* 


(group) 


T(C)T I T(w)":" (C) (C)+ "}" 


(min gap) 


(integer) 


(max gap) 


(integer) 


(window size) 


(integer) 



Table 1 - Definitions of a Pattern 



[0044] For example, a pattern {A}1 - 3{B} is a pattern with two groups 
15 each having a single category A and category B, respectively. The minimum gap 
and maximum gap between the two groups are 1 and 3 time units, respectively. 
In a pattern {2 : A, B}6 - 8{3 : B, A, B}13 - 15{C}, there are three groups {2 : A, B}, 
{3 : B, A, B}, and {C}. The first group has a window size of 2 and contains two 
event categories, one each of event category A and B. The second group has 
20 one category A and two copies of category B. The third group only has a single 
event category, category C. 

[0045] From this, a general pattern may be described by Formula [1]: 

{ Wi:E 1f1 E 1in i} Ch-pi{ u) 2 : E 2.1.---. E 2 ,n2}a 2 - - P 9 -i{ oo g : E gf1 ,..., E g , ng } [1] 

25 

[0046] An occurrence of this pattern in a given sequence is a subset of 
event instances in the sequence {ei,i,..., ei.ni, e^i,..., e 2 , n2 ,..., e g ,i,..., e g , ng }. As 
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such, a pattern occurrence is a set of event instances satisfying particular 
structural and temporal constraints such as described in formula [2] to [5]. 



C(e ifj ) = E itj ; when 1 < i < g and 1 < j <n g . [2] 

rnaxj T(ei,j) - minjT(ej,j) < wf, when 1 < j < g. [3] 

maxj T(e i+1J ) - minjT(ejj) < Pi; when 1 < i < (g -1). [4] 

minj T(e i+ ij) - maXjT(eij) > ar, when 1 < i < (g -1 ). [5] 



[0047] The event categories in the pattern and the categories are in 
one-to-one correspondence. Event instances that are mapped to the same group 
occur in the corresponding window size. The time gap between any two event 
instances mapped to two consecutive groups is less than or equal to the 
specified maximum gap and no less than the specified minimum gap. For each 
event instance in the occurrence, an event instance matches an event category 
in the pattern when the event instance is mapped to the category in the one-to- 
one correspondence. 

[0048] Two different occurrences of the same pattern are disjoint when 
the intersection of the two sets is empty. A set of occurrences of a pattern is a 
disjoint occurrence set when any two occurrences in the set are disjoint. 

[0049] As an example, in sequence 100 of Fig. 1 , a pattern in the 
illustrated temporal sequence may be identified. Counting the number of the 
occurrences may be based on a rule definition, as sequence 100 includes 
overlapping occurrences of some patterns. For example, one overlapping pattern 
is an event A followed by an event B in one to three time units. This pattern is 
defined as {A}1 - 3{B}. The occurrences of A may be identified and then a 
determination of whether an instance of B follows within the constrained temporal 
region of 1 to 3 time units. However, as there are multiple instances of B 
satisfying the pattern constraint, the method defines sets instances of "disjoint 
occurrences" to include all legitimate occurrences including those with 
overlapping events or categories. However, different occurrences in the same 
set do not share event instances. 
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[0050] A sequence may contain many different sets of disjoint occur- 
rences. A maximum cardinality c of these sets of disjoint occurrences is the 
number of such disjoint occurrences and may be referred to as the occurrence- 
based frequency or o-frequency. The maximum cardinality c of all disjoint 
5 occurrence sets of a pattern in a sequence is defined and identified as a number 
or count of pattern occurrences in the sequence. The maximum cardinality c is a 
function of the sequence and the pattern and does not require an additional 
parameter such as a sliding window size. The maximum cardinality c ensures 
that patterns occur more often when a temporal constraint is relaxed. The 

10 maximum cardinality c is an occurrence-based frequency of a pattern in a data 
sequence for all disjoint occurrence sets in that sequence. If a sequence does 
not contain at least one occurrence of a pattern, the maximum cardinality c is 
zero. If a sequence contains at least one occurrence of the pattern, there is a 
disjoint occurrence set with the maximum cardinality c. As such, the maximum 

15 cardinality c has a lower-bound of zero and an upper-bound equal to the number 
of occurrences of the event category in the pattern having the least number of 
occurrences. As will be discussed further below, the maximum cardinality c may 
also provide for the estimation of the expected number of pattern occurrences or 
the expected maximum cardinality. This estimation of the expected number of 

20 pattern occurrences may also be compared to the counted occurrences to 
identify a surprising pattern. 

[0051] The method provides, in one implementation, for determination 
of the maximum cardinality c by determining the quantity of the temporal patterns 
in sequences such that each event instance is used no more than once in 

25 counting pattern occurrences. For example, sequence 100 has six non-empty 
sets of disjoint occurrences of pattern {A}1 - 3{B}. These six non-empty sets are 
identified in Table 2 along with a count of the disjoint occurrences. 

30 
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1 


{{A@11,B@14}} 


1 


2 


{{A@13, B@14}} 


1 


3 


{{A@13, B@15}} 


1 


4 


{{A@13, B@16}} 


1 


5 


{{A@11,B@14}, 


2 




{A@13, B@15}} 




6 


{{A@11,B@14}, 


2 




A@13, B@16}} 





5 

Table 2 - Occurrences and Count of Disjoint Occurrences 

[0052] As shown in Table 2, the maximum cardinality c of pattern {A}1 - 
3{B} in sequence 100 is two. As shown in Fig. 1 and Table 2, in occurrence set 

10 number 5, event category A instance 106 at time period 1 1 is only counted once. 
Event category A instance 108 at time period 13 is only counted once in 
determining the maximum cardinality. Similarly, event category B instance 1 10 at 
time period 14 is only counted once and event category B instance 1 14 at time 
period 15 is only counted once. As such, occurrence set number 5 has two 

1 5 disjoint occurrences, and the maximum cardinality is at least two. As shown in 
Fig. 1 , event category A instance 1 1 6 at time period 1 6 is not followed by any 
other event category B instance, and there are only two other event category A 
instances 106 and 108. As such, the maximum cardinality is at most two. 

[0053] As the method defines maximum cardinality, all occurrences of a 

20 pattern P are occurrences of a relaxed pattern P r of pattern P. Pattern P may be 
relaxed by several methods. A relaxed pattern P r may have one or more 
categories dropped or removed without dropping all categories from any group. 
Also, a relaxed pattern P r may have the first or the last group dropped. In the 
alternative, a relaxed pattern P r may include an increased window size for one or 

25 more groups. Similarly, a relaxed pattern P r may include increasing one or more 
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maximum group gaps. In yet another relaxed pattern P r , one or more minimum 
group gaps may be decreased. 

[0054] As such, the maximum cardinality of all disjoint occurrence sets 
of the relaxed pattern P r is greater than or equal to the maximum cardinality of 
5 the original pattern P. For a given sequence the maximum cardinality of a pattern 
P with additional or narrower constraints, maximum cardinality is less than or 
equal to maximum cardinality of the broader pattern. For example, for pattern 
{A}1 - 3{B} in sequence 100, a more constrained pattern {A}2 - 3{B} also has a 
maximum cardinality of two or less. 

10 [0055] In operation, one or more implementations provide for a method 

of determining the maximum cardinality c of a pattern P in a sequence s. In a 
method 200 of Fig. 2, one implementation of the method receives an input 
sequence s in operation 202 and a defined pattern P in operation 204. In 
operation 206, the method identifies all simple occurrences of the pattern within 

15 the subsets of event instances in sequence s. One or more disjoint occurrence 
sets are identified from the set of simple occurrences in operation 208. The sets 
of identified disjoint occurrence sets are counted in operation 210. The maximum 
cardinality c of pattern P in sequence s is output in operation 212. 

[0056] This method may be applied easily to short data sequences. 

20 However, for large amounts of data and/or long sequences, such an approach 
has an exponential computing cost 

[0057] In other embodiments and implementations, a system and 
method determines distinct occurrences of a pattern in one or more sequences of 
time-stamped event instances. The method includes determining a maximum 

25 cardinality of disjoint occurrences of the pattern in the one or more sequences. 
The method may be embodied in computer executable instructions within a 
computer readable medium. The system or computer executable instructions 
include storing the sequence and means for defining the pattern. Also included is 
determination of a maximum cardinality of disjoint occurrences of the pattern in 

30 the sequence. 

[0058] In one implementation, the method includes identifying a single 
disjoint occurrence set that has maximum cardinality. The method includes 
scanning the temporal data sequence along the time line. This scanning may be 
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in the forward direction as illustrated, or may be in the reverse direction. The 
method may be a sequential flow with various process loops. In the method, the 
loop index V is a pointer indicating the current event under evaluation with a 
process operation of one or more of the process loops. Index i is therefore 
5 between zero and the total number of events "m" in the sequence. The index i is 
moved back and forth during the method. The method transitions when the index 
i is greater than the total number of events m in the sequence. The number of 
matched occurrences "c" is updated during this process based on the evaluation 
and determinations of the method operations. 

1 0 [0059] Referring to Fig. 3, a method 300 illustrates another 

implementation of a method for determining the maximum cardinality of pattern P 
of operation 304 in sequence s of operation 302. Sequence s from operation 
302 and pattern P from operation 304 are input into operation 306 to initialize 
method 300. The method continues in operation 308 with the determination of 

15 the occurrences of the pattern P in sequence s. This determination may be 
performed by mapping pattern recursively to the events and categories of 
sequence s or may be performed by using a pattern template or other method, 
one or more of which are described herein. Event instances in the determined 
occurrences are identified in operation 310 and the disjoint occurrences are 

20 identified using the determined occurrences from operation 308 and the identified 
event instances in the occurrences from operation 312. A disjoint occurrence 
count is provided from operation 312 to operation 318. 

[0060] In operation 314, instances contained in the disjoint occurrence 
set of operation 312 are removed from the sequence under consideration. 

25 Operation 31 2 ensures that further analysis of the sequence does not utilize 
event instances or event category instances that have already been utilized in 
counting a disjoint occurrence set. In an alternative implementation, operation 
312 could flag the event instances in lieu of removing them. Operation 316 
provides a looping analysis function by comparing the current loop counter index i 

30 to the total number of events m in the sequence. When the index i is less than or 
equal to the total number of events m, operation 316 loop method back to 
operation 308 for further determination and analysis. Once index I is greater than 
the total number of events m, operation 316 breaks the looping and provides an 
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indication to operation 318 that the process is complete. Operation 318 sums the 
total number of disjoint occurrences and provides operation 320 with the count. 
Operation 318 may also provide a listing of each disjoint occurrence set to 
operation 320 or as another output of method 300. Operation 320 provides the 
5 maximum cardinality of pattern p in sequence s as an output of method 300. 

[0061] Referring now to Fig. 4, method 400 is one implementation of 
the method of determining the occurrences of a pattern P in a sequence. As 
noted above, pattern P is defined as having one or more groups. This is one 
implementation of the method operation 308 of Fig. 3. Method 400 receives the 

10 initialization from operation 306 that includes the one or more sequences to be 
analyzed and the pattern P. In operation 404, event of index i is matched to a 
group j where j is an index number for the group under consideration. 

[0062] Operation 406 increases loop index i by one and operation 408 
checks to see if group j is matched. If group j is not matched, operation 408 

1 5 routes method 400 back to operation 404 for further analysis. If group j is 
matched in operation 408, operation 410 checks to determine if the group is 
within the window of group j. If it is not, operation 410 routes the method to 
operation 412 where the loop index i is set to h, a secondary index. Secondary 
index h may be specified to be less than index i thereby moving the method 

20 backward from the previous index i. Next operation 414 removes all matches in 
group j and returns the method to operation 404 for further analysis. 

[0063] If in operation 41 0 group j is matched and within the window of 
group j, the method continues in operation 416 where group j and group j minus 
one are determined to be within the maximum time gap defined by pattern P. If 

25 the time gap between group j and group j minus one is not within the maximum 
time gap of pattern P, operation 416 routes the method to operation 418 where 
the index i is reset to the smallest index. In operation 420, group index j is 
decreased by 1 so that in the next looping operation group j minus 1 is reworked 
or analyzed. 

30 [0064] If in operation 41 6 the time gap between group j and group j 

minus one is within the maximum time gap for pattern P, the method continues to 
operation 422 where the occurrence count is incremented. The incremented 
occurrence count and/or the occurrence from operation 422 are provided to 
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operation 310 that provides for the identification of event instances in the 
identified occurrences as described above with regard to operation 310 in Fig. 3. 
Additionally, operation 422 routes method 400 to operation 424 where the current 
loop index i is compared to the total number of events m in sequence s. If index i 
is less than or equal to the total number of events m, the method is looped back 
to operation 404 for further analysis. If however, index is greater than the total 
number of events m, the method is routed to operation 426 where the 
occurrences and occurrence count are output in operation 426. 

[0065] Referring now to Fig. 5, method 500 illustrates another 
implementation for determining the maximum cardinality of a pattern in a 
sequence. The maximum cardinality of pattern P is determined in a given 
sequence s:{ei, e 2 , . . . , e m }, where events e m are temporally ordered by 
timestamps. Method 500 starts with an initialization of variables for the number 
of matched occurrences c, the number of matched groups j, and the index i of the 
event or instance in sequence s in operation 502. 

[0066] To find an occurrence of a pattern, the method matches on a 
group-by-group basis, To match a group, multiple events may be required as 
may be specified by the pattern definition. Events within the sequence are 
matched to the events within a group of the pattern P in operation 504. 
Operation 506 checks whether group j is fully matched and whether index i is 
within the maximum events m in sequence s. Event index V is verified to be less 
than or equal to an "m" total number of events in sequence "s" such that the 
index i is within the current sequence s under consideration. As will be 
discussed, this verification is associated with the looping within operation 504 and 
operations 506 to 516. In these method operations, looped analysis of the 
sequence is performed and index i is updated as a function of one or more loops 
of the method of operation 506. 

[0067] If group j is not fully matched and index i is within the total 
number of events m in the sequence in operation 506, the current event i is 
matched against group j in operation 508. In operation 510, index i is 
incremented to the next event. As illustrated here, index i is increased by 1 so 
that the next event is matched in a forward analysis loop. However, in an 
alternative implementation, index i may be initiated to total number of events m in 
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operation 502 and decreased by 1 in operation 510 thereby providing for a 
backwards analysis loop. 

[0068] All the events that result in a match group are determined and 
identified in operation 508. Once all categories in a group are matched to some 
5 event instance, the window size constraint on the group is checked in operation 
512. Operation 512 determines whether all event categories in group j are 
matched and analyzes whether the window width u)j for group j is within range or 
whether it is violated. If the constraint is violated, the method reverses or moves 
backward by a unit specified by the window size coj. The method attempts to map 

1 0 the same group to try to identify another match. 

[0069] If not within the group window width constraint, the matched 
events in the group j are discarded in operation 516. Similarly, if the group j is 
not fully matched and the window size ooj is not violated in operation 512, the 
method loops back to operation 506 for further matching. The match process is 

1 5 looped back to operation 506 at a determined index later than the previous 

matching attempt. The methods within operation 504 are repeated until index i is 
out of range, i.e., index i is greater than the total number of events m in sequence 
s, and group j cannot be matched, or a match of group j is identified within the 
group window constraint. 

20 [0070] When operation 516 is complete, no event category in group j is 

set as matched. Additionally, a new starting point for index i is established in 
operation 514 as secondary index h to rematch group j. The starting point h will 
be set greater than index i because the window width ujj for group j was violated. 
The method operation 514 provides that the method does not repeatedly find the 

25 same set of matching events that violate the window width u)j constraint. 

Operation 512 checks the matched events for compliance with the window width 
constraint as defined by the group definition. 

[0071] The method illustrated by sub-operations 506 to 516 within 
operation 504 continuously loop until either a group j is matched or the index i is 

30 greater than the number of events m in sequence s as determined in operation 
506. When group j is fully matched or index i is greater than total number of 
events m, the method checks to determine if group j cannot be matched and 
whether index i is greater than the total number of events m in operation 518. 

20 
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[0072] The number of matched occurrences is returned at operation 
520 when group j cannot be matched and index i is greater than the total number 
of events m in sequence s. As such, operation 520 reports the number of disjoint 
occurrences of pattern P in sequence s, e.g., the maximum cardinality c. 

[0073] if in operation 518, group j can be matched and loop index i is 
less than or equal to the total number of events m in sequence s, the method 
continues at operation 522. That is, in operation 518 if there is at least one 
match of group j, then the process goes to operation 518 for further matching. In 
operation 522, the method checks whether the most recent group j is within the 
time-constraint with respect to the previous group, group j minus one, e.g., 
whether the maximum time gap p ri is violated. When the most recent group j 
violates the maximum time gap (3 n , the previous group j minus one is re-matched 
in addition to the current group j. 

[0074] If the maximum time gap (3 n is violated, the starting position for 
matching the next group is determined in operations 524 and 526. In operation 
524, the number of matched groups is decreased by 1 , and the previous group 
and any groups following the previous group are re-matched. If the group is 
outside of the range, e.g., it is too far away from its preceding group per operation 
522, the match for this group and the match for the preceding group are 
discarded. In this case, the method moves backward by a unit of the maximum 
gap, and attempts to match the preceding group as in operation 524 and 526. 

[0075] When the maximum time gap (3 ri is not violated in operation 
522, the method continues to operation 528. When the time constraints of 
window width and time gaps are satisfied, operation 528 considers the partial 
match of the pattern from the first group to the current group to be valid. When 
the group is within range from the preceding group and the group is successfully 
matched, the method moves forward to skip the minimum gap in operation 528, 
and to skip the first item matched in this group in the previous matched pattern 
occurrence per operation 528, to match the next group. 

[0076] In operation 532, the index i is advance by the minimum time 
gap a] from the latest time in group j. In optional operations 530 and 532, the 
index i is advanced so that index i is greater than the first event of group j in the 
previously matched occurrence of pattern P. Optional operation 530 provides, in 
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one implementation, for optimization of the method for some patterns within some 
sequences. In operation 532, index i is increased to a time of an event of group j 
that is greater than the earliest time in group j of the last match occurrence. In 
operation 534, the number of matched group j is increased by 1 and the process 
5 continues to operation 536 for matching of pattern P as a function of the matched 
group j. In addition, if in operation 530 there are no disjoint occurrences or fully 
matched patterns, then the method goes to operation 534 where the step index i 
is increased by one. 

[0077] In operation 536, the method matches the groups to determine 

10 the occurrence of the whole pattern P such that the entire pattern definition of 
groups and window gaps are matched, e.g., no window width or group gap 
violations. If there are no violations, then an occurrence of pattern P is 
determined. In operation 538, the number of occurrences c is increased by one 
and the matched events either are removed from the sequence or are flagged in 

15 operation 540. The group number index j is reset to zero in operation 542 and 
may be looped back to operation 504 and specifically to operation 506 for further 
analysis. As an optional optimization method, once a pattern occurrence is 
identified in operation 536, all matched instances are removed from the sequence 
in operation 540. Event instances occurring before the first one in the matched 

20 pattern occurrence are ignored in operation 544 and the method loops back to 
operation 506 to find the next occurrence. This looping continues until group j is 
not fully matched and index i is greater then the m number of events in sequence 
s. 

[0078] The method of removing or flagging of the matched events in 
25 operation 540 within the matched occurrence of pattern P ensures that a second 
or future matching does not include the same previously matched events. As 
previously discussed, the method of Fig. 5 provides for the determination of the 
maximum cardinality that is the number of discrete occurrences of pattern P in 
the one or more sequences. As such, no two discrete occurrences may have 
30 one or more events in common. However, two matched discrete events may 
temporally overlap. As such, removal or flagging of the events within the 
matched occurrence in operation 540 provides for further matching of events, 
groups, and patterns that do not re-use a previously matched event within the 
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temporal sequence but enables further matching of temporally overlapping 
occurrences. 

[0079] After operation 540, the method loops back to identify another 
occurrence of pattern P. Operation 536 checks to identify whether the identified 
5 groups comply with the pattern definition. 

[0080] At the conclusion of operation 536, when at least one 
occurrence is identified, the matched group counter j is reset to zero in operation 
542 and index i is reset for another loop. Index i is reset in operation 544 to an 
event immediately after the earliest event that was previously matched and 
10 removed or flagged. This optional method provides for an efficient method by 
looping back to a previous match rather than starting from the beginning. 

[0081] Operation 504 is looped until loop index i is greater than the total 
number of events m in sequence s. Once index i is greater than the maximum 
number of events m, operation 504 directs the method to operation 518 that then 
1 5 directs the method to operation 520. Operation 520 provides for the output of the 
then current number of occurrences, e.g., maximum cardinality c. Outputting in 
operation 520 may include reporting, storing, transmitting, etc. the maximum 
cardinality c or o-frequency of pattern P in sequence s or in a plurality of 
sequences s. 

20 [0082] Operation 504 addresses the matching of a group. The method 

does not start at the same time index i to match the same group more than once 
due to the methods of operations 514, 524, and 532. As such, the method 
provides for a maximum number of loops of O(mg) in operation 504. For each 
loop in operation 504, the group match operation 506 has a method cost of 

25 O(mn), where n is the number of categories in the pattern. Operations 518, 522- 
544 have a method cost O(m). Hence the overall method cost is less than or 
equal to 0(m 2 gn). 

[0083] Optionally, further operations not illustrated may be provided, 
some of which provide for optimization in particular situations. For example, if 
30 there are event categories in the pattern with very few occurrences in the 

sequence, such occurrences may be searched first to find pattern occurrences 
around such occurrences without requiring the matching of another part of the 
pattern. 
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[0084] For example, in an alternative implementation of operation 516, 
part of the matched occurrences satisfying the window width ojj constraints may 
be reused rather than resetting the matching process for the j-th group. Other 
optional operations may be added that may provide for optimization of the 
5 method when addressing particular patterns and/or addressing particular 
sequences of temporal events. 

[0085] The maximum cardinality provides for the determination or 
identification of a pattern within time-stamped sequence that includes overlapping 
temporally-related events. Such a method and system provides, in an example 

10 application, for the identification of related maintenance events. Once identified, 
improved maintenance practices may be prepared that provide for reduced 
equipment maintenance costs and improved equipment reliability. 

[0086] Application of some implementations of the method described 
herein provide for the identification of new patterns rather than a simple search or 

15 retrieval and counting of an existing or known pattern. 

[0087] Some implementations and embodiments may address the 
recurrence of a pattern in the same sequence of temporally-related events. 
Counting the occurrence frequency of a pattern in a sequence by determining the 
number of discrete recurrences of each pattern within the sequence. 

20 [0088] While other data mining methods are very costly to enumerate 

all sets of occurrences for each pattern against each sequence in a database of 
time stamped events, embodiments of the present disclosure provide for 
reductions in the computational costs for data mining patterns in complex time- 
stamped event sequences. 

25 [0089] In another implementation, the method provides for an 

estimation of the maximum cardinality. The method applies a probability 
assumption to determine the expected quantity of disjoint occurrences of the 
pattern as a function of various characteristics and/or parameters. These may 
include the total quantity for each event instance or category of event instances 

30 within a sequence and as included in the pattern. Also, these may include the 
maximum and minimum time gaps as defined by the patter and the maximum 
time length of the sequence. 
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[0090] In the exemplary implementations and embodiments discussed 
herein, for each type of event in a sequence, the method and system assumes 
that all instances are uniformly distributed for the time range of the sequence. 
Such example method also assumes that all events and all types of events are 
5 independently distributed. However, this is for illustration purposes as the 

method contemplates other event distributions may also be assumed and utilized 
in a similar manner. 

[0091] The method assumes that the time of each instance in a 
sequences is a random variable, independent from other instances and uniformly 
10 distributed over the sequence time range. The method specifies that X x , y is a 
random variable of the number of occurrences of pattern {A}a - p{S} having x 
occurrences of A and y occurrences of B and that where y > x > 0. The method 
assumes that Sx, y; b is the probability that X x , y = k, i.e., ^ X) y;b := Pr(X x , y = b). The 
method includes p y := Si, y: i and q y := ^ 1>y;0 . As such, p y = (1 - p) y and q y = 1 - p y = 
15 (1-pp. 

[0092] For any give x, y, and b, where 0 < b < x: 

^x,y;0 = ^ 1,y;0 S x-1,y;0 = q y £ x-1,y;0 = q X y [6] 

£x,y;b = £ 1,y;0 £x-1,y;b + ^1,y,1 £x-1,y-l,b-1 = Py 5x-1.y-1;b-1 \7] 

20 £x,y;x = % 1,y;1 ^ x -1,y-1;x-1 = Py €x-1,y-1;x-1 [8] 

Mi.y = Sl.y.1 = Py [9] 

Mx.y = Exp(X x ,) = p y + Py|Jx-1,y-1 + q y JJx-1,y [10] 

[0093] For any given event or category A, the probability p that any 
25 category or event B associated with category or event A is a discreet occurrence 
of pattern P where the probability p is equal to or greater than zero and less than 
or equal to one. In the method, the probability estimate p x>y of the value of the 
expected maximum cardinality is determined. Probability estimate p x , y may be 
determined by determining the mean p 1>y for y instances of event B. In one 
30 implementation, the probability estimate p XiV is determined by a recursive 

method. For example, to determine p 3 ,4, both p 2 ,3 and p 2 ,4 must be determined. 
Similarly, to determine p 2 ,3, both p 1t2 and p 13 must be determined. However, in 
practice such a recursive method is very costly. 

25 
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[0094] In one alternative implementation, the method estimates the 
mean p Xiy as a function of a bound on the mean of the expected maximum 
cardinality. The method assumes there are x unique instances of pattern P that 
each consisting of q number of event B's. The method also assumes there are y 
5 instances of event B. The y instances of event B is greater than or equal to the 
product of the q sum of the random variables and the x number of incidences of 
event A, e.g., y number of instances of event B is greater than or equal to the 
produce of q and the x number of instances of event A. 

[0095] In such an implementation, the method determines two 

1 0 boundary values of the mean of the expected maximum cardinality. To obtain an 
upper bound when counting the pattern {A}a - p{8}, all event B's within pattern P 
are reused. To obtain the lower bound, one implementation of the method uses 
the determined number b of pattern P and separately y - qb number of event B's. 
The determined b number of pattern P and y - qb number of event B are used to 

15 count the combination of pattern P and event B. As such, in one implementation 
where an event B follows pattern P, the method determines a mean of the 
expected maximum cardinality to be bound by Mx, y and max b=1 . x (p kiy . nb ). 

[0096] In an alternative implementation, the bounds of the mean of the 
expected maximum cardinality are determined as a function of an estimation 

20 precision objective p. The estimation precision objective p is an integer or 
fraction of the real mean |j r such that the estimated bounds are within the 
estimation precision objective p, i.e., the upper bound |j + minus the lower bound 
p" should be less than or equal to the estimation precision objective p. As such, 
where two times the loop counter a plus one is less than x number of 

25 occurrences of a category in sequence s, then the upper bound p + can be refined 
to meet the objective. Where two times the loop counter a is less than x, then the 
lower bound p~ can be refined to meet the objective. 

[0097] The expected maximum cardinality may be determined by 
determining the lower bound p'and upper bound p + to the mean jj of the 

30 expected maximum cardinality. One such method 600 is illustrated in the 
flowchart of Fig. 6. As with the above, the pattern P is {A}a - p{B} is only 
described here for illustration purposes only. It should be noted, however, similar 
operations may be applied to other patterns. 
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[0098] Method 600 determines the lower bound p~ and the upper bound 
p + for the mean of the expected maximum cardinality as a function of an 
estimation precision objective p of the real mean p r . When the x number of 
instances of A is small, the maximum number a of method loops is also small, 
5 and the bound may not be tight or within estimation precision objective p of the 
real mean p r - However, when index i is larger than a predetermined threshold 
imax, the estimation yields very tight bounds within the estimation precision 
objective p of the real mean p r . Such a threshold i max may be a function of 
various parameters including the length of the sequence, the categories or events 

10 in the sequence, and the pattern P time gaps. 

[0099] For this method, the x number of occurrences of event A, the y 
number of occurrences of event B, the probability p (that derives q from one 
minus the probability p), and estimation precision objective p are provided as 
inputs into the method and the lower bound p" and the upper bound p + are 

1 5 provided as outputs. 

[00100] In the implementation illustrated in Fig. 6, method 600 
determines an estimate for the lower bound p~ and the upper bound p + for the 
expected maximum cardinality c. In operation 602, the method determines an 
estimation coefficient cp and an alternating adjustment factor A. In operation 604, 

20 a base estimation of the mean of the expected maximum cardinality 5 is 

determined. In operation 606, an incremental estimation parameter qj e for the 
event or category of events is determined and an incremental estimation 
parameter qj Q is set to a first instance of the incremental estimation parameter qj^ 
[00101] Next, operation 608 determines an estimation adjustment 

25 parameter y. Loop counter a is set to one in operation 610. In operation 612, an 
initial lower bound p~ is determined as a function of adding the incremental 
estimation parameter ip e to the base estimation 5 of the mean of the expected 
maximum cardinality. An initial upper bound p + is determined by adding the null 
incremental estimation parameter k\) 0 and the estimation adjustment parameter y 

30 to the base estimation 5 of the mean of the expected maximum cardinality. 

[00102] After the initial upper bound p + and lower bound p" are 
determined in operation 612, the method refines the upper bound p + or lower 
bound p* by a looping analysis. Operation 614 checks the upper and lower 
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bounds to determine if their combination is within a predetermined range and 
variance. Operation 614 determines whether the relative difference between the 
upper and lower bounds (as defined by p + minus p') is smaller than a predefined 
estimation precision value p around the real mean Mr- If the difference is less 
5 than or equal to the product of the estimated precision value p and the real mean 
p r , no further refinement of the upper bound p + or lower bound p~ are required. As 
such, the looping breaks to route the method to operation 616 and the initial 
values of the upper bound p + and lower bound p" are provided as the range of the 
expected maximum cardinality. If however, the upper and lower bounds is 

10 greater than the product of the estimated precision value p and the real mean p r , 
method 600 refines the lower bound p" and/or the upper bound p + in a looping 
process until neither can be further refined. 

[00103] For example, in one implementation, further refinement to the 
upper bound p + is determined to not be necessary when two times the a number 

1 5 of loops is equal to or greater than x number of occurrences of a category in the 
sequence. Similarly, the lower bound p" may be evaluated to determine whether 
further refinement is desired. This occurs when the sum of two times the number 
of loops a and one is equal to or greater than x number of occurrences of a 
category in the sequence. 

20 [00104] In order to refine the upper and lower bounds, estimation 

coefficient <p k(2 as a function of the current loop k and the estimation index z are 
determined in operation 618. Additionally, alternating adjustment factor A k for the 
current loop k is determined in operation 618 for two values, the first at two times 
the current loop index a and the second at the sum of two times the current loop 

25 index a and 1. 

[00105] The method checks the current loop index a in operation 620. If 
the current loop index a is a value such that the product of two and the current 
loop index a is greater than or equal to x number of occurrences of a category in 
the sequence, the loop breaks. The method routes to operation 616 where the 

30 current values of the upper and lower bounds are reported as the range of the 
expected maximum cardinality. 

[00106] If however, the product of two and the current loop index a is 
less than x number of occurrences of a category in the sequence, the refined 
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factors of operation 618 are used in operations 622 and 624 to refine the lower 
bound p\ The incremental estimation parameter nj e for event e is determined in 
operation 622 and the lower bound pis re-determined as a function of the re- 
determined incremental estimation parameter qj e in operation 624. 
5 [00107] Operation 626 re-checks the variance condition against a 

predefined variance threshold or the estimation precision objective p. Operation 
626 checks the range of the upper and lower bounds to determine if it is within a 
predetermined range and variance. Operation 626 determines whether the 
relative difference between the upper and lower bounds (as defined by p + minus 

10 p~) is smaller than a product of the predefined estimation precision value p and 
the lower bound p~. If the difference is less than or equal to the product of the 
estimated precision value p and the lower bound p~, further refinement of the 
upper bound p + or lower bound p~ is not required. As such, the looping breaks to 
route the method to operation 616 and the initial values of the upper bound p + 

15 and lower bound p" are provided as the range of the expected maximum 

cardinality. If however, the upper and lower bounds is greater than the product 
of the estimated precision value p and the lower bound p", method 600 refines the 
upper bound p + . 

[00108] Operation 628 checks the current loop index a. If the current 
20 loop index a is a value. The sum of one plus the product of two and the current 
loop index o is compared to the x number of occurrences of a category in the 
sequence. If the sum is greater than or equal to the x number of occurrences, 
method 600 breaks and routes to operation 616 for reporting of the current values 
of the upper and lower bounds as the range of the expected maximum 
25 cardinality. 

[00109] If however, the sum is less than the x number of occurrences, 
operation 630 updates the incremental estimation parameter ip G - The upper 
bound p + is updated in operation 632 as a function of the updated incremental 
estimation parameter hj q of operation 630. The method continues to the next loop 
30 by indexing the loop index by one in operation 634. After operation 634, the 
method loops back to operation 614 to determine whether further refinement is 
desired by checking whether the re-determined bound range is within the 
estimation precision range. 
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[00110] Generally, the method illustrated in Fig. 6 determines the 
estimated upper bound p + and lower bound p~ in a three process method. First, 
initialization derives an initial lower bound p~ and upper bound p + through 
operation 608. A second process refines the lower bound p~ and than the upper 
5 bound p + until neither of them can be further refined by operations 620 and 628, 
or the bounds are tight such that their relative difference is equal to or smaller 
than the estimation precision objective p. Third, the refined estimated upper 
bound p + and lower bound p~ are provided as outputs in operation 616. It should 
be understood, however, that in another implementation, the order of refining 

10 upper bound p + and lower bound p~ may be reversed, or both may be refined in 
the same method operation. 

[001 11] Another implementation of the method according to the present 
disclosure is disclosed in an algorithm form in Appendix B. The method provides 
for a determination of the lower bound p'and upper bound p + ofthe mean of the 

15 expected maximum cardinality of pattern P in sequence s. The minimum bound 
p and maximum bound p + are first estimated and then refined to within a 
predetermined estimation precision objective. The method with the estimation 
precision objective a ensures a bounding of the mean of the expected maximum 
cardinality utilizing an incremental and iterative looping process. 

20 [00112] In another implementation, a method for determining the 

estimated bounds of expected maximum cardinality is a function of estimation 
and bound tightness relationships as provided in formula [1 1] to [25]: 
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Ai£ (P/, z q k(2+1) [17] 



5 = x-(qy- x+1 (1-q x ))/(1-q) [18] 
= 0 [19] 
Yi = q 2y - 2x+3 (1-q) [20] 



2a- 



MJ2a + 1 = ( d 2a,1 + d 2 o + 1,l) [21] 

2a-2 

Y2a = (dza-U+dsa.l) [22] 



/=i 



Yk = X (A,£ cpu(q (z+1)(k+1) -q x(z+1) )/(1-q 2+1 ) [23] 

'=1 2=0 

1 5 |J + = 5 + HJ2G-M + Y2a+1 [24] 

M" = 6 + l|l2o+A 2 o [25] 

20 [001 1 3] Some implementations of the method apply one or more of the 

relationships defined in formula [1 1] to [25] to determine the lower bound p'and 
upper bound jj + of the mean of the expected maximum cardinality. One such 
implementation is illustrated in method 700 of Fig. 7 and another implementation 
is illustrated in Appendix C. As with the method of Fig. 6, method 700 illustrates 

25 one implementation of the method of the present disclosure by addressing 
exemplary pattern P is {A}a - (3{B}. 

[001 14] Method 700 determines the lower bound p" and the upper bound 
p + for the mean of the expected maximum cardinality as a function of an 
estimation precision objective p of the real mean p r . In operation 702, the method 

30 determines an estimation coefficient cp and an alternating adjustment factor A. In 
operation 704, a base estimation 6 of the mean of the expected maximum 
cardinality is determined. Base estimation 5 may be determined by formula [18] 
as indicated in operation 704. In this operation, base estimation 6 is a function of 
the predetermined probability p of the occurrence, the x number of occurrences 

35 of event or category A, and the y number of occurrences of event or category B. 
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Of course if patter P defined others events or categories, their counted 
occurrences would also be a factor in this determination. 

[00115] In operation 706, an incremental estimation parameter i|j e for 
the event or category of events is determined and an incremental estimation 
5 parameter ip 0 is set to a first instance of the incremental estimation parameter 
Incremental estimation parameter ijj e may be defined per as indicated in 
operation 606 to be equal to zero. Incremental estimation parameter qj Q is set to 
be equal to incremental estimation parameter ijj-i. Incremental estimation 
parameter 14^ is defined in one implementation by formula [20] or in operation 706 

10 as being a function of the predetermined probability p of the occurrence, the x 
number of occurrences of event or category A, and the y number of occurrences 
of event or category B as indicated. 

[00116] An estimation adjustment parameter y is determined in 
operation 708 or in another implementation by application of formula [25]. 

1 5 Estimation adjustment parameter y is a function of alternating adjustment factor 
A, the predetermined probability p of the occurrence, and the x number of 
occurrences of event or category A. 

[00117] The loop counter a is set to one in operation 710. An initial 
lower bound p~ is determined in operation 712 as a function of adding the 

20 incremental estimation parameter qj e to the base estimation 5 of the mean of the 
expected maximum cardinality. In addition, an initial upper bound p+ is 
determined in operation 712 by adding the odd incremental estimation parameter 
ip 0 and the estimation adjustment parameter y to the base estimation 5 of the 
mean of the expected maximum cardinality. 

25 [00118] After the initial upper bound p + and lower bound p~ are 

determined in operation 712, the method checks to determine if further 
refinement of the upper bound p + and/or lower bound p* is desirable. The 
method checks the range of the upper and lower bounds about the real mean p r 
to determine if the range is within a predetermined estimation precision range. 

30 The predetermined bound range of the difference between the upper and lower 
bound (as defined by p + minus p") is compared in operation 714 to the estimation 
precision range (as defined as the product of the predefined estimation precision 
value p and the real mean p r ). If the bound range is less than estimation 
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precision range, than method 700 determines that the initial upper bound |j + and 
lower bound p~ are sufficiently tight about the real mean and further refinement of 
the upper bound p + and/or lower bound p~ is not required. As such, method 700 
breaks and routes to operation 716 where the current values of the upper bound 
5 p + and lower bound p" are provided as the range for the expected maximum 
cardinality. 

[00119] If however, operation 714 determines that the bound range is 
greater than or equal to the estimation precision range, method 700 refines the 
upper bound p + and the lower bound p~ in operations 718 to 734 in a looping 

1 0 process. The looping process of operations 71 8 to 734 continue until the bound 
range as determined in operation 714 is less than the estimation precision range. 

[00120] For example in one implementation, further refinement to the 
lower bound p~ is determined to not be necessary when two times the a number 
of loops is equal to or greater than x number of occurrences of a category in the 

1 5 sequence. Similarly, the upper bound p + is tested to determine whether further 
refinement is desired. This occurs when two times a number of loops of the 
method plus one is equal to or greater than x number of occurrences of a 
category in the sequence. When either the upper bound p + or lower bound p~ are 
determined to not require further refinement, the looping operation of method 700 

20 stops and the then current values of p + and p'are reported in operation 716 as 
the maximum and minimum bounds of the mean of the expected maximum 
cardinality. 

[00121] In order to refine the upper and lower bounds, estimation 
coefficient cp k>z as a function of the current loop k and the estimation index z are 
25 determined in operation 718. Additionally, alternating adjustment factor A k for the 
current loop k is determined in operation 718 for two values, the first at two times 
the current loop index a and the second at the sum of two times the current loop 
index a and 1. 

[00122] The method checks the current loop index a in operation 720. If 
30 the current loop index a is a value such that the product of two and the current 
loop index a is greater than or equal to x number of occurrences of a category in 
the sequence, the method breaks. The method routes to operation 716 where 
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the current values of the upper and lower bounds are reported as the range of the 
expected maximum cardinality. 

[00123] If however, the product of two and the current loop index a is 
less than x number of occurrences of a category in the sequence, the refined 
5 factors of operation 718 are used in operations 722 and 724 to refine the lower 
bound The incremental estimation parameter ip e for event e is determined in 
operation 722 and the lower bound p~ is re-determined as a function of the re- 
determined incremental estimation parameter iy e in operation 724. 

[00124] Operation 726 re-checks the variance condition against a 

10 predefined variance threshold or the estimation precision objective p. Operation 
726 checks the range of the upper and lower bounds to determine if it is within a 
predetermined range and variance. Operation 726 determines whether the 
relative difference between the upper and lower bounds (as defined by p + minus 
p') is smaller than a product of the predefined estimation precision value p and 

1 5 the lower bound p\ If the difference is less than or equal to the product of the 
estimated precision value p and the lower bound p", further refinement of the 
upper bound p + or lower bound p~ is not required. As such, the looping breaks to 
route the method to operation 716 and the initial values of the upper bound p + 
and lower bound p" are provided as the range of the expected maximum 

20 cardinality. If however, the upper and lower bounds is greater than the product 
of the estimated precision value p and the lower bound p", method 600 refines the 
upper bound p + . 

[00125] Operation 728 checks the current loop index a. If the current 
loop index o is a value. The sum of one plus the product of two and the current 

25 loop index a is compared to the x number of occurrences of a category in the 
sequence. If the sum is greater than or equal to the x number of occurrences, 
the looping operation of method 700 breaks and routes to operation 716 for 
reporting of the current values of the upper and lower bounds as the range of the 
expected maximum cardinality. 

30 [00126] If however, the sum is less than the x number of occurrences, 

operation 630 updates the incremental estimation parameter ijj 0 . The upper 
bound p + is updated in operation 732 as a function of the updated incremental 
estimation parameter hj q of operation 730. The method continues to the next loop 
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by indexing the loop index by one in operation 734. After operation 734, the 
method loops back to operation 714 until one of the predefined non-refinement 
criteria are met or exceeded. 

[00127] Generally, the method illustrated in Fig. 7 determines the 
5 estimated upper bound p + and lower bound p~ in a three process method. First, 
initialization derives an initial lower bound p~ and upper bound p + through 
operation 708. A second process refines the lower bound p~ and upper bound p + 
until neither of them can be further refined by operations 720 and 728, or the 
bounds are tight such that their relative difference is equal to or smaller than the 
10 estimation precision objective p. Third, the refined estimated upper bound p + and 
lower bound p" are provided as outputs in operation 716. It should also be 
understood that, in another implementation, the order of refining upper bound p + 
and lower bound p" may be reversed, or both may be refined in the same method 
operation. 

1 5 [00128] As described, method 700 assumes that both event category A 

and B are independently and uniformly distributed. However, in other 
implementations, the method may assume and utilize events distributed under 
other distributions such as a Poisson distribution. In such cases, the methods 
disclosed herein for estimation of the minimum and maximum bounds of the 

20 mean of the expected maximum cardinality may be adopted for use with Poisson 
distributions or any other definable probability distribution. 

[00129] The methods of Fig. 5, Fig. 6, Fig. 7, Appendix B, and Appendix 
C describe exemplary implementations for the method of estimating the expected 
maximum cardinality of a pattern within in sequence. Such methods for 

25 determining the mean of the expected maximum cardinality generally provide for 
improved operational performance of data mining systems and reduced 
computational costs. 

[00130] Such an estimated maximum cardinality is useful in data mining 
for identifying a surprising or interesting pattern within one or more sequences. 

30 The identified surprising pattern may be a pattern that was not previously known 
or expected and therefore may be an indication of a change to the relationship of 
events, or may be the identification of a new pattern or relationship between 
events. 
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[00131] A surprise pattern in a sequence may be identified by 
determining an estimated of the expected maximum cardinality and comparing it 
to the determined maximum cardinality of the pattern in the one or more 
sequences. Such a comparison is a measure of the dependency of the 
5 correlated events in a sequence. Such a measure may be referred to as lift. 

[00132] A method 800 as illustrated in Fig. 8 is one implementation for 
identifying a surprise pattern. One or more sequences may be input into method 
800 in operation 302. A pattern P is input into method 800 in operation 304. A 
maximum cardinality of pattern P in sequence s is determined in operation 802. 

10 Operation 802 operates one or more of the methods described above for the 
determination of the maximum cardinality. In Operation 804, the expected 
maximum cardinality is estimated. The expected maximum cardinality may be 
determined by one or more of the methods described above such as the 
estimation of the upper and lower bounds to the mean of the expected maximum 

15 cardinality. The determined maximum cardinality from operation 802 and the 
estimated expected maximum cardinality from operation 804 are provided to 
operation 806. Operation 806 identifies a surprise pattern by application of one 
or more determinations. Such determinations may include a comparison, a 
trending, an analysis, or otherwise. 

20 [00133] For example, in one implementation, where the determined 

maximum cardinality of a pattern in one or more sequence deviates significantly 
from the expected value, a surprising pattern may be identified. Operation 806 
may identify a surprising pattern when there is a large maximum cardinality and 
where the maximum cardinality differs from an expected value by more than a 

25 threshold level. In an exemplary implementation, a large maximum cardinality 
may be a maximum cardinality of greater than about 10. The threshold level for 
the difference may be about 20 percent. In such a case, operation 806 may 
identify a pattern, as a surprise pattern. 

[00134] Operation 806 may also identify a surprising pattern when the 

30 determined or counted maximum cardinality is small and the absolute difference 
between the determined maximum cardinality and the expected value is greater 
than a small maximum cardinality threshold amount. A maximum cardinality may 
be small if, for example, it is less than or equal to about 10. A small maximum 
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cardinality threshold amount may be a threshold of greater than about 30 
percent. 

[00135] In another implementation, method 900 as illustrated by the flow 
chart in Fig. 9 also provides for the identification of a surprise pattern. One or 
5 more sequences may be input into method 900 in operation 302. A pattern P is 
input into method 900 in operation 304. A maximum cardinality of pattern P in 
sequence s is determined in operation 802. Operation 802 operates one or more 
of the methods described above for the determination of the maximum 
cardinality. Operation 902 receives the sequence from operation 302, the pattern 

10 from operation 304, and the determined maximum cardinality from operation 802 
to estimate the expected maximum cardinality. The expected maximum 
cardinality may be determined by one or more of the methods described above 
such as the estimation of the upper and lower bounds to the mean of the 
expected maximum cardinality. The determined maximum cardinality from 

15 operation 802 and the estimated expected maximum cardinality from operation 
902 are provided to operation 806. Operation 806 identifies a surprise pattern by 
application of one or more determinations as described above. 

[00136] As discussed, the determination of the maximum cardinality of a 
pattern in one or more sequences does not require extra parameters such as a 

20 sliding window size and does not depend on a particular algorithm or set of 

algorithms. The method is monotonic as patterns do not occur less frequently 
when their constraints are relaxed. The method provides for the determination of 
the maximum cardinality of a pattern within one or more sequences that may be 
implemented in a system using an efficient greedy algorithm. Additionally, the 

25 method estimates an expected maximum cardinality that can be compared to the 
determined maximum cardinality. From this, a pattern in a sequence may be 
evaluated to identify a surprising pattern. 

[00137] The method estimates the maximum cardinality under different 
assumptions of data sequences thereby providing for identifying surprising 

30 patterns. The expected maximum cardinality is determined under an inde- 
pendent uniform distribution assumption. The independent uniform distribution 
assumption may reflect knowledge with regard to the dependency between 
events or the determined maximum cardinality. In practice, other pattern or 
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distribution assumptions may be applied using the method. For example, a 
Poisson distribution may be adapted to the method based on the domain 
knowledge, and the knowledge of the event timing distributions within sequence 
s. 

5 [00138] Persons skilled in the art will understand that the method 

disclose herein may be implemented in hardware or software and may be defined 
in software code, as a flow or decision chart, as one or more forms of a computer 
program product, and/or in an algorithm, such as a simple greedy algorithm. One 
or more computing systems that include programming means and/or computer 

10 readable medium with computer instructions for operating in accordance with the 
methods and operations of the methods described herein are included with the 
scope of the present disclosure. 

[00139] Fig. 10 illustrates one embodiment of a high level block diagram 
for a system for performing one or more embodiments of the present disclosure. 

15 A data mining computer system 1000 includes one or more processors (not 
shown) that may be configured with computer instructions for data access and 
data mining of a sequence containing time-stamped events or categories and 
analysis according to the present disclosure as described above. Data mining 
computer 1000 may be configured to determine the maximum cardinality of a 

20 pattern in a sequence. Data mining computer system 1000 may also determine a 
mean of the expected maximum cardinality of a pattern in a sequence and 
compare that with the determined maximum cardinality to identify a surprise 
pattern. 

[00140] Data mining computer 1000 may include one or more databases 
25 1002. Database 1002 includes one or more sequences of time-stamped events 
and/or event categories. Database 1002 may be any type of database including 
a simple compilation of data in a simple spreadsheet or word processing file. 
Data mining computer 1000 may include a memory 1004 that is any type of 
memory used in a computing or information processing system. A data access 
30 program, utility, or module 1006 accesses data stored in the data files of 
database 1002 and/or in memory 1004 for analysis. 

[00141] A user input device 1008 receives user input associated with the 
selection of the sequence or sequences to be analyzed. Additionally, user input 
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1008 may receive a pattern definition to be data mined including one or more 
parameters or characteristics of the pattern. These may include one or more 
events, categories, time gaps, and or windows. User input 1008 may be directly 
associated with data mining computer 1000 or may be remote whereby user 
5 definable criteria and/or variables are provided via a communications link or 
channel (not shown). User input 1008 may be any type of input for example, 
another computer system, a personal data assistant, a pointing device, a 
keyboard, a memory, etc. 

[00142] Data access module 1006 may provide a sequence and a 

10 pattern to a maximum cardinality module 1010 that includes a data mining sub- 
module 1012 and a disjoint occurrence counting sub-module 1014. Maximum 
cardinality module mines the sequence data to identify disjoint occurrences of the 
target pattern and counts the disjoint occurrences to provide a maximum 
cardinality of the pattern in the sequence to an output module 1022. 

15 [00143] Data access module 1006 may provide the sequence data and 

the pattern to estimation and surprise pattern identification module 1016. 
Estimation sub-module 1018 receives the sequence, the pattern, and may 
receive the determined maximum cardinality from counting sub-module 1014. 
Estimation sub-module 1018 determines the expected maximum cardinality 

20 based on a probability distribution assumptions as discussed above. In some 
embodiments, estimation sub-module 1018 may determine bound p + and lower 
bound m~ of the mean of the expected maximum cardinality. 

[00144] A comparison sub-module 1020 receives the estimated 
maximum cardinality from the estimation sub-module 1018 and the determined 

25 maximum cardinality from counting sub-module 1014. Comparison sub-module 
1018 performs a comparison of the estimated maximum cardinality and the 
determined maximum cardinality. Comparison sub-module 1018 may include 
one or more methods, criteria, and/or thresholds to provide for analysis of the 
pattern as a function of the estimated maximum cardinality and the determined 

30 maximum cardinality. Comparison sub-module 1018 may identify a surprising 
pattern as discussed above and provide its data and results to output 1022. 
Additionally, comparison sub-module may collect surprising or interesting 
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patterns, rank them according to one or more criteria, filter them based on one or 
more filter characteristics, and/or analyze the patterns further. 

[00145] Data mining computer output 1022 provides the received data 
from counting sub-module 1014 and/or comparison sub-module 1020 to storage 
5 in a memory 1024, to a display 1026, to a communication interface or facility 
1028 for transmission to a remote system, to a local or remote printer 1030, 
and/or to a report module or generator 1032. Such outputs 1022 may include a 
visual or graphic representation of one or more surprising patterns, their rank, or 
various categories of patterns based, at least in part, on filtering one or more 

1 0 using a filter characteristic. 

[00146] Generally, data mining computer 1000 may be configured with 
computer executable instructions on a computer readable medium for performing 
one or more of the methods disclosed above. 

[00147] In operation, some embodiments of the present disclosure 

1 5 provide for reduced computational processing requirements and cost of data 

mining a pattern from one or more sequences. For example, the inventor tested 
a data set of 25 sequences each containing 105 events. Each of the 105 events 
were independently, uniformly and randomly assigned to one of 10 categories 
and distributed on a time line. Several tests were preformed to determine the 

20 impact on the computational cost for data mining and to determine the 
robustness of the method. 

[00148] One embodiment of the present disclosure was tested to 
analyze the data set and to count different patterns in the data set. The method 
was implemented in Java and the test was run within Sun JRE 1.4.1 on a laptop 

25 with a 1 .2GHz CPU and 51 2MB of main memory. All 25 sequences were loaded 
in main memory and the method in the form of software code was initiated. 

[00149] The average runtime and count was tested for patterns with 
different temporal constraints. Two types of patterns were considered: {C1}1- 
p{C2}1- (3{C3} and {oj: C1, C2, C3}1- P{C4}, where C1, C2, C3, C4 are the first 

30 four event categories, and co = p varies from 10 to 106. In other words, the 

expected number of events in the search region (to match the next event in the 
same group or the next group) varies from 0.1 to 10,000. It was determined that 
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patterns with multiple categories in a group are more difficult to count and more 
sensitive to the increase of the search region. 

[00150] The method was found to scale sub-linearly with the search 
region size. The runtimes ranged from 6.4 milliseconds to 18.8 milliseconds for 
the first pattern and 9.6 milliseconds to 709.4 milliseconds for the second pattern. 
This reflected a considerable improvement over other methods. 

[00151] Another test determined the average runtime and count for 
patterns with different lengths. In this case, two types of patterns were 
considered, {C1}1- p{Ck} where each group contains only a single category, and 
{u> : C 1 , . . . , Ck-1}1- p{Ck}. C1,..., Ck are the first k event categories, p = 1000, 
i.e., the search region (for a category) contains approximately 10 events in the 
sequence. In this case, the test demonstrated that the method provides for a 
linear scaling as a function of the length of patterns. 

[001 52] The scalability as a function of the sequence size of one 
embodiment of the method was also tested. In this test, three additional data 
sets were generated using similar settings except that the total number of evens 
(103, 104, and 106 respectively) and time range ([1, 105], [1, 106] and [1, 108] 
respectively) in each sequence are different for different data sets. Again, time 
period [t + 1 , t + 100] was expected to contain a single event on average for all 
sequences. Two patterns were tested for averaged runtime and count: {C1}1- 
p{C2}1 - p{C3} and {u> : C1, C2„ C3,}1 - p{C4} where C1,..., C4 are the first four 
event categories and to = 0 = 1000. The total runtime for counting a single 
pattern in 25 sequences in the smallest data set was less than 1 millisecond. In 
testing, both the count and the runtime increased linearly as a function of the 
sequence size. Such as linear increase was an improvement over other methods 
that demonstrated a quadratic increase. 

[00153] In this case, for pattern {C1}1- (3{C2}1 - P{C3}, the time required 
for processing was 0, 1.2, 11.6, and 114.2 milliseconds for sequences sizes 103, 
104, 105, and 106, respectively. The determined counts were 28.6, 301.4, 
2,980.3, and 29,817.8, respectively. For pattern {co : C1, C2„ C3,}1 - p{C4}, the 
time required for processing was 0, 1.6, 16.8 and 167.8 milliseconds for 
sequences sizes 103, 104, 105, and 106, respectively. The determined counts 
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were 18.4, 186.2, 1858.4, and 18,696.0, respectively. These results represent a 
significant reduction in required processing time. 

[00154] By application of the methods disclosed herein to other data 
mining methods and systems, similar improvements in data mining efficiency 
5 may be provided. 

[00155] While the examples and descriptions are generally described 
with regard to a single sequence, this is only exemplary and is not intended to be 
limiting. It should be understood by one skilled in the art, the disclosed method 
and system may also determine maximum cardinality of a one or more patterns in 

1 0 two or more sequences. For example, a single sequence may be the 

maintenance events and records for one aircraft. However, two or more of the 
sequences of aircraft records may represent the maintenance records of a fleet of 
aircraft. As such, the present disclosure may identify a pattern that is common 
across two or more aircraft or across the fleet of aircraft or may identify a 

1 5 surprising pattern of maintenance related events across the fleet. 

[00156] It is further to be understood that the methods or operations 
described herein are not to be construed as necessarily requiring their 
performance in the particular order discussed or illustrated. It is also to be 
understood that additional or alternative operations may be employed. 

20 [00157] When introducing aspects of the present disclosure or 

embodiments thereof, the articles "a", "an", "the", and "said" are intended to 
mean that there are one or more of the elements. The terms "comprising", 
"including", and "having" are intended to be inclusive and mean that there may be 
additional elements other than the listed elements. 

25 [00158] While various embodiments have been described in whole or in 

part, those skilled in the art will recognize modifications or variations that might 
be made without departing from the inventive concept. The examples illustrated 
in the present disclosure and are not intended to limit it. Therefore, the 
description and claims should be interpreted liberally with only such limitation as 

30 is necessary in view of the pertinent prior art. 



42 



Boeing Ref. 02-1381 
HDP Ref. 7784-000707 



10 



15 



Substitute Specification (Clean Version) 



APPENDIX B 



1 Method initialization 

1.1 Determine (pi, z and ^ 

1.2 Determines 

1.3 Set hj e = 0, ip 0 = ifi 

1 .4 Determine Yi 

1.5 Assume a = 1 

1 .6 Assume u" = 5 + ip e and u+ = 5 + iy 0 + Yi 

2 if u + - u" £ pu", then break, if not then 

2.1 Determine q> kz and A k for k = 2o, 2o + 1 

2.2 If 2o > x, then go to 3 

2.3 Update iy e 

2.4 Update u" 

2.5 If p + - u " < pm", then go to 3 

2.6 If 2o + 1 > x, then go to 3 

2.7 Update ip 0 

2.8 Update u + 

2.9 Set a = a + 1 

2.10 Go to 2 
20 3 Output u + and u" 
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Substitute Specification (Clean Version) 

Boeing Ref. 02-1381 
HDP Ref. 7784-000707 

APPENDIX C 

1 Method initialization 

1.1 Determine (p 1j2 and A.,: 

<Pi,o = 1 

<Pi.i = - 1 
A, = q 2 <*- x+1 > 

1.2 Determined: 5= x - (q y x+1 (1-q x ))/(1-q) 

1.3 Set ip e = 0, ip 0 = : where 
MJ1 = q 2y - 2x+3 (1 - q) 

1 .4 Determine Yi 

Y^Ai ((q 2 -q x )/(1-q)-(q 4 -q 2x )d-q 2 )) 

1.5 Set a = 1 

1 .6 Set p" = 5 + ijJe and |j+ = 5 + ip 0 + yi 

2 while (p + - p" ^ pp~), then break, if not then 

2.1 Determine (p kiZ and A k for k = 2a, 2a + 1: 

<Pk,o = £ ' 2 

z=l 

9k.z = -((p k -i,z-i)/(1-q z ) 
A k = -A^d-qV^ 1 

2.2 If (2a > x), then go to 3 

2.3 Update ip e : 

increase^ by 2 £ [A z £ (p z .w(q (2o+1)(w+1) +q 2o(w+1) )] 

z=l w=0 

2.4 Update p~: Set p~ = 6 + qj 2o + Y20 

2.5 If (p + - p- < pu"), then go to 3 

2.6 If (2a + 1 > x), then go to 3 

2.7 Update ip 0: 

increase^ by £ [A z £ cp z , w (q (2a+1) < w+1 >+q 2o < w+1 >)] 

2.8 Update p + : Set p + = 5 + ip 2o +i + Y2a+i 
2.8 Set a = a + 1 

3 Output p + and p" 
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